RADIO SCIENCE Journal of Research NBS/USNC-URSI 
Vol. 69D, No. 1, January 1965 



Alteration of the Electron Density of the Lower Ionosphere 
With Ground-Based Transmitters 



Pietro P. Lombardini 

Moore School of Electrical Engineering, University of Pennsylvania. Philadelphia, Pa. 

(Received May 15, 1964; revised July 28, 1964) 

Using computational techniques as described by Carleton and Megill, and Megill to 
calculate the difference between rate of attachment of electrons to 2 and the ionization rate 
in air as a function of electric field, the power required to create breakdown in the ionosphere 
has been calculated. These calculations are compared with extrapolated microwave break- 
down data and found to agree satisfactorily. 

The propagation of the primary beam is altered by the electrons created by it. r l 'his 
gives origin to a nonlinear problem, which is very difficult to deal with in the general case. 
However, a steady state approximate solution has been found, considering diffusion in one 
dimension. The solution is applied to the design of a possible atmospheric experiment. 



1. Introduction 

The possibility of altering the characteristics of 
the lower ionosphere utilizing ground-based trans- 
mitters has been considered by several investigators. 
In particular, studies have been made concerning 
the achievement of a decrease in electron density 
[Molmud, 1964J, an increase in the electron collision 
irequency [Farley, 1963] and the excitation of air 
glow [Bailey, 1959 and Megill, 1964]. Little atten- 
tion, however, has been given to utilizing ground- 
based transmitters to achieve an increase of electron 
density. The aim of this paper is to investigate the 
parameters involved in obtaining such an increase. 

The first step of the present work involves the 
definition of breakdown as created by a high-frequency 
electromagnetic wave in a boundless, slightly ionized 
homogeneous medium. Experimental data 1'or these 
conditions of breakdown are not available; however, 
the electric field necessary to initiate breakdown 
(breakdown threshold) may be obtained in either 
of the following two ways. One way consists in 
performing direct laboratory measurements of break- 
down in finite enclosures, and then scaling the results 
to the open air case. The other way consists of 
computing the electron distribution function in the 
mixture of gases considered (N 2 , 2 , and O) and in 
the presence of the electric field. This distribution 
function computed by utilizing the experimental data 
for the various cross sections involved is in turn 
used to evalute the electron attachment and 
ionization frequencies (i.e., v a and v t respectively). 
The breakdown threshold, E b (defined as that field 
for whicn v L =v a ) may therefore be calculated. 

The breakdown threshold is shown to decrease 
for certain altitudes at a faster rate than that 
specified by the inverse of the distance law which 



applies to VHF waves in the far field of the trans- 
mitter. From this, it is concluded that breakdown 
can be achieved if the proper power is available. 
Furthermore, it is shown that once breakdown is 
obtained, the steady state electron density is critically 
dependent upon the difference of the altitude rate of 
changes of the breakdown field and the inverse of the 
distance field. 

Because of this dependence on the altitude rate of 
change of the exciting field, greater electron densities 
may be achieved using a sufficiently large antenna 
such that breakdown occurs in the near field. 

A characterization of the attachment-ionization 
process is shown to indicate an altitude for which 
the optimum breakdown efficiency exists. 

2. Characterization of the Breakdown 
Threshold Electric Field 

A VHF radio wave traveling upward in the D or 
in the lower E region is assumed to be idealized as a 
plane wave propagating in a weakly ionized, quasi- 
homogeneous plasma. The breakdown conditions 
applicable may therefore be identified with those 
pertaining to an electric field of constant orientation 
and amplitude in a weakly ionized air. In general, 
the breakdown condition of an electric discharge is 
defined by the electric field (breakdown threshold, 
E h ) for which equilibrium is achieved between the 
gain of electrons due to ionizing collisions, and the 
loss due to the various mechanisms which remove 
electrons from the discharge region (e.g., attachment, 
recombination and diffusion). Because of the ab- 
sence of walls, and the weak ionization for the case 
considered, the breakdown is predominantly con- 
trolled by attachment, while diffusion has a negligible 
role. Conditions of this kind are rarely encountered 
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in the laboratory where the walls confining the 
discharge play a considerable role in removing 
electrons. Brown and Rose [1957] have shown 
that the experimental data of microwave breakdown 
threshold in air can be characterized by a curve of 

#eff • 

— — versus the product pL, where E ett is the " effec- 
tive" breakdown electric field defined by the 
formula, 



or the rate of decrease of electrons due to attachment, 



E ett =E t 



V^ + co 2 



(1) 



where. 



v= average electron collision frequency 
co = radian frequency of wave 
#=air pressure 
L= breakdown gap spacing. 

This curve flattens out for increasing values of pL 
approaching the asymptotic value of 3X10 3 V/m 
x (mm of Hg). Introducing this value into (1), 
considering a collision frequency as a function 
of pressure given by 5.3X10 9 p(c/s) [MacDonald, 
Gaskell, and Gitterman, 1963], one obtains the 
following simple equation relating the rms value of 
the breakdown threshold field in air with the pressure, 



E b 



=3X10 8 pJi+(- 2 Y- 

^\ ^V5.3X10V 



(2) 



Specification of breakdown conditions in open space 
may also be obtained utilizing a different approach. 
That is, the energy distribution function of the 
electrons under the influence of an oscillating electric 
field in a mixture of N 2 , 2 , and O may be calculated. 
In this approach, the experimental data introduced 
are the collision cross sections of the electrons with 
the above mentioned components of air. Carleton 
and Megill [1962], have given a detailed analysis of 
the problem of determining the electron energy 
distribution function under these conditions. Megill 
has constructed a computer program for the IBM 
704 of Raytheon, which gives the electron energy 
distribution function, f(u), for the interval of 
electron energy, u, from 0.01 to 20.0 eV, and having 
neutral particle, density, temperature, geomagnetic 
field, frequency, and the electric field intensity as 
parameters. The solution applies to a situation 
in which the total number of electrons does not 
change. Providing that the rate of change of the 
total number of electrons is small, and incorporating 
perturbation techniques, the solution j(u) may be 
used for evaluating the rate of increase of electrons 
due to ionizing collisions. 



Vi = ~ 
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a (u)du 



(j) m f™ 'v l y(u)dv 



(3) 
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(4) 



where N iy <j U) and a iA are the concentration, ioniza- 
tion cross section and attachment cross section for 
the ith component of air considered (i.e., N 2 , 2 , 
and O) and m is the electron mass. In the actual 
computations, the only attachment phenomenon con- 
sidered is the dissociative attachment, 2 +e-^0+0~. 
Ionization and attachment collision frequency ob- 
tained in this way may be characterized as a function 
of the electric field intensity E, holding all other 
parameters constant. v a grows rapidly for moderate 
fields and quickly tends to saturate, while v t starts 
a quick rise at higher field intensities. The crossing 
of the two curves (vi=v a ) corresponds to the break- 
down threshold E b . Using the atmospheric model 
presented in the Handbook for Astronautical Engi- 
neering [Koelle, 1961] extensive computations have 
been made for the following choice of parameters: 
/=50 Mc/s, geomagnetic fields 5 X10 _5 Wb/m 2 , and 
concentrations and temperatures corresponding to 
the height interval, 55 to 105 km. The corre- 
sponding values of E b as a function of altitude are 
shown in figure 1 by the solid curve. In the same 
figure, the dotted curve shows the E b altitude varia- 
tion as obtained using (2). The good agreement 
between the two methods of approach is apparent 
from this figure. The collision frequencies given 
by (3) and (4) are required for evaluating the rate 
of change of the total number of electrons given by, 



dN , 



■Va)N. 



(5) 



Computed values of vt— v a are plotted in figure 2 
as a function of the excess of electric field, AE, 
and the same range of parameters previously speci- 
fied for figure 1. The relationship between v x —v a 

and -^- is approximately linear. 

In figure 3, v i -—v a is plotted versus height for 
constant electric field above threshold breakdown. 
This figure shows that there exists an altitude at 
which the creation of new electrons occurs with a 
maximum efficiency. This optimum height cor- 
responds to the condition for which the average 
collision frequency approximately equals the fre- 
quency of the electric field; i.e., 



v{h, E b )~f. 



(6) 



Figure 4 shows the relationship between frequency 
and altitude as obtained by (6), using for v, Mac- 
Donald, Gaskell, and Gitterman [1963] approxima- 
tion. To each of these frequencies corresponds a 
breakdown threshold electric field, which may be 
determined using (2). 
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Figure 1. Breakdown threshold electric field versus altitude for a frequency of 50 Mc/s. 
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Figure 2. TTie difference of the ionization and attachment collision frequencies as a function of the percent increase of excess electric 

field above breakdown. 
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Figure 3. The difference of the ionization and attachment collision frequencies as a function of altitude for the field strengths. 

indicated. 



3. Maximum Attainable Increase of Electron 
Density- 
In the previous section, the conditions required for 
obtaining breakdown have been examined. The 
problem is now that of determining the maximum 
electron density which may be achieved using radio 
waves. For this purpose, consider a VHF wave 
propagating vertically upward. Under normal iono- 
spheric conditions, a wave in this band of frequencies 
is negligibly absorbed by the ionosphere. However, 
when the electric field strength of the wave becomes 
larger than the breakdown threshold, the electron 
density will increase and hence the wave may suffer 
an exponential attenuation. In free space, the 



electric field of the wave may be expressed as: 
E(x)=E $(x) 



(7) 



where E is the amplitude of field intensity at an 
altitude, h , from the ground, where the source is 
assumed to be located. Here, x is the local vertical 
coordinate in a system whose origin is at a height 
h , and $(x) represents the functional variation in 
the absence of absorption. In the "far field" of an 
antenna, 



$(x)=- 



1 



1 + 



K 



(8) 
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Figuee 4. Frequency of maximum ionization efficiency versus altitude. 
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Along the axis of a spherical antenna excited in 
phase and focused at A , *(x) is given by, 



sin 



$(x)=- 



\8\hf; 






(9) 



where D is the antenna diameter and X the wave- 
length. 

If the propagation takes place in ionized air, the 
wave is subjected to exponential absorption. When 
the electron density is sufficiently low, the differential 
power absorbed along the wave path is, 



im=-»<*> 



2 (eE(x)) 2 
3 



m 



/» 00 

U 3/2 a(u) ^C v w/ du 



&/(«) 



r« 1/! /'(«)rfo 



(10) 



where N(x) is the electron density at the altitude 




A + x and a (u) is given by, 



x(u)-- 



v 2 +u z 



(11) 



disregarding the earth's magnetic field. Proceeding 
to integrate (10) in the interval x=0 to the arbitrary 
distance x, 

E(x)=E $(x) exp f~- f X A(E, x)N(x)dx\ 

(12) 

where the quantity A (E, x) is a slowly varying 
function of E and x. 

In order to simplify the problem, let us for the 
moment disregard diffusion. Assume that for the 
radio waves propagating upward with some spatial 
variation, <t> (x) (e.g., inverse of the distance law) 
breakdown conditions are reached at the altitude hi. 
The physical situation may be explained from 
figure 5 where the spatial variation of the electric 
field of the wave and the breakdown threshold 



ALTITUDE- 



Figuke 5. Sketch of the variation of the free space field strength (dashed line), and of the breakdown threshold (solid line) versus 

altitude. 
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variation are plotted as a function of altitude. 
These two curves cross at two altitudes, hi and h 2 . 
Consider a point located at the altitude hi-\-x, where 
the spatial variation field curve lies above the 
breakdown curve. Here the intensity of the wave 
will change with time because of the attenuation 
due to the time increasing- electron density. Hence, 
assume the spatial variation curve of figure 5 repre- 
sents the established situation at time, t=Q. It 
follows at this point (hi+x) that the electron density 
will start increasing exponentially with time accord- 
ing to the solution of (5), 



where 



N(x,t)=N exp[(v i -v a )t] 



(13) 



where iV is the previously existing electron density. 
It has already been pointed out that the quantity 
(vf—Va) is critically dependent upon the quantity 
AE=E—E b (see fig. 2). Substituting (13) into (12), 



AE=E $(x) exp 

- £ A(x)N exp [(vi-VaWdxl-Etix). (14) 

Equation (14) does not take into account any 
mechanism limiting the growth of electron density 
as a function of time, such as losses due to diffusion 
or recombination. Therefore, for increasingly large 
times, satisfaction of (14) requires that AE ap- 
proaches zero. In the limit condition, 



Lim P A(E, x)N(x, t)dx= -In J^- (15) 
Nmax=Lim N(x, t) 



<o0(z) 

' d4>(x) dE b (x) ^~ 
dx — dx 



(16) 



An important observation of (16) is that the 
steady state distribution of N is critically dependent 
upon the slopes of the field and the breakdown 
threshold curves. 



4. Diffusion Problem 

Equation (16) was obtained ignoring diffusion. 
The rigorous solution taking diffusion into account 
is difficult. However, for steady state, if the electric 
field and the density of the neutral components of 
air are slowly varying functions of the coordinates 
x, the electron balance may be approximated by 
the equation: 



D(E, x) 



d 2 N{x)_ 



D(E, x)= 



r v?' 2 v- i (u)j(u)du 

_ Jo 

3 f °° 

- m u l,2 j(u)du 



(is) 



dx 2 



-N(x)(v t -v a ) 



(17) 



is the diffusion coefficient which is assumed to be 
a slowly varying function of the altitude and the 
electric field. 

Equations (12) and (17) describe fairly well the 
steady state situation for N (x) and E (x). The 
system can be solved using the boundary conditions 
suggested by physical considerations (see fig 6). 
The boundary conditions will thus be chosen at two 
altitudes hi and h 2 , one below and the other above 
the crossing of the breakdown threshold. For this 

altitude, the value of N and -7— are assumed to 

be known. Examining (12) and (17) we note that 
while A (E, x) and D (E, x) are slowly varying 
functions of E and x, the function p t — v a may undergo 
large changes for small changes of either variable 
(as shown in fig. 6). Thus, (17) is very sensitive to 
small changes of E. This renders the system of (12) 
and (17) highly nonlinear and difficult to deal with. 
This system has been solved numerically by Alber- 
toni, Bocchieri, and Daneri [1963] with the help of an 
IBM 7090 computer. 

5. Numerical Examples 

Two cases are considered here which employ the 
results of the previous section. The first case deals 
with the "far field' ' of the transmitter. The param- 
eters selected are, 

/=50 Mc/s 

hi=71 km 

# & =196V/m 

In figure 7 is plotted the electron density distribution 
as obtained using (13) (curve on the right). Also 
plotted (curves on the left; note the change in scale) 
are the solutions of the system of (12) and (17). 
The two curves, which correspond to the different 
"trial" parameters chosen in the numerical integra- 
tion, represent the physical solution as long as they 
coincide. Beyond the point of bifurcation they fail 
to indicate a physical solution. The region of coin- 
cidence could be prolonged with a better choice of 
parameters. However, the data in figure 8 shows 
that beyond the region considered, the approximation 
obtained neglecting diffusion (13) is adequate. 

The second case refers to a hypothetical trans- 
mitting antenna large enough such that the break- 
down field is located at its focus. Here it is assumed 
that the antenna is a spherical basin whose surface 
is radiating in phase. The law of variation of the 
field intensity near the focus is thus assumed to be 
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Figure 6. Variation af the difference of the ionization and attachment collision frequencies as a function of the field strength. 
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Figure 7. Electron density distribution for case 1 [f = 50 Mc/s E b = 196 V/m, far field of transmitter]. 



given by (9) where, 
ttD 2 



= 1.7XHT 2 /m. 



8\Ai 2 

The other parameters are chosen such that in absence 
of electronic absorption, the electric field radiated 
reaches a maximum of 200 V/m at the altitude of 78 
km, where the focus of the antenna is assumed to be 
located. In figure 8 is plotted the electron density 
variation with altitude obtained using (16). Also 
plotted are the solution of the system of (12) and (17) 
which is valid as a physical solution up to the point of 
bifurcation. 

6. Power Requirements 

No mention has been made yet concerning the 
power required for achieving breakdown. The 
ground-based upward radiating transmitter needed 
for obtaining the previously mentioned breakdown 
fields must have a power-gain product PG which for 
free wave propagation is given by: 



PG=A7rh 2 E b 2 / Vo 



(19) 



where 770 is the wave impedance in space. 

The power P fed into the antenna is determined 
once the gain G of the antenna is known. For the 
case in which the ionization takes place in the far 
field of the transmitting antenna, the antenna dimen- 
sions are limited by the condition : 



h>2D 2 /\. 



(20) 



Assuming a paraboloid antenna of diameter D, the 
maximum antenna gain can not exceed the limiting 
value : 

£ m ax = 7T 2 V2A. ( 21 ) 

On the other hand, when the ionization takes place 
in the near field of the antenna, and the antenna is 
focused, the value of the gain is only limited by 
practical factors like construction accuracy and cost. 

In table 1 are tabulated breakdown transmitter 
parameters for the conditions of altitude and fre- 
quency given by figure 4. For the far field case the 
maximum gain in the table is obtained using formula 
(21). For the near field case, the assumed gain at 
focus in the table is based on gains achieved by far 
field antennas for the same frequency. 
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Figure 8. Electron density distribution for case 2 [f = 50 Mc/s, E(hi)= 200 V/M, E b = 156 V/m, near field of transmitter). 



Table 1 





Frequency 




Far field case 


Near field case 


Break- 


for 
optimum 


Power-gain 
product 














down 


ionization 


PG 


Maximum 




Assumed 




altitude 


efficiency 




gam 


Power 


gain 
at focus 


Power 


km 


Mc/s 


Watts 




Megawatt 




Megawatt 


60 


1, 000 


1.2 1015 


1.0 106 


1,200 


7. 5 106 


160 


65 
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3. 8 1014 


5.5 105 
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3. 7 106 
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70 
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1.0 101^ 
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1. 9 106 


53 


75 
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2. 10'3 


1.3 105 
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7. 5 105 


18 
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66 


3. 105 


12 


Note that the powers listed above are relative to the condition needed for start- 


ing breakdown. Therefore, in order to obtain an increase in electron density, 


the power required should be larger than the ones tabulated according to the 


square of 
ithe same 

i 


the ratio of 
iltitude. 


the field use 


d to the field 


needed for i 


tarting bres 


ikdown at 
J 



7. Conclusions 

It has been shown that a radio wave radiated by 
a ground-based transmitter may achieve breakdown 
in the D or lower E ionospheric regions. For each 
altitude, a frequency exists at which there corre- 
sponds an optimum ionization efficiency. The maxi- 
mum electron density which can be obtained criti- 
cally depends upon the altitude rate of change of both 
the breakdown field and the free space transmitted 
field strength. This characteristic automatically 
limits the maximum obtainable electron density 
which may be achieved in the far field of an antenna. 
This limitation may, however, be overcome by 
having breakdown take place in the near field. 
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